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Magnetic Catalysis and Oscillating Effects in the Nambu — Jona-Lasinio Model 

at Nonzero Chemical Potential [] 

K.G. Klimenko Q 

Institute for High Energy Physics, 142284 Protvino, Moscow region, Russia 

Phase structure of the four dimensional Nambu - Jona-Lasinio model has been investi- 
gated in two cases: 1) in nonsimply connected space-time of the form R 3 x S* 1 (space coordi- 
nate is compactified and the length of the circle S 1 is L) with nonzero chemical potential fi 
and 2) in Minkowski space-time at nonzero values of fi, H, where H is the external magnetic 
field. In both cases on phase portraits of the model there are infinitly many massless chirally 
symmetric phases as well as massive ones with spontaneously broken chiral invariance. Such 
phase structure leads unavoidably to oscillations of some physical parameters at L — > oo 
or H — > 0, including magnetization, pressure and particle density of the system as well as 
quark condensate and critical curve of chiral phase transitions. Phase transitions of 1st and 
2nd orders and several tricritical points have been shown to exist on phase diagrams of the 
model. 



1. Introduction § 

The concept of dynamical chiral symmetry breaking (DCSB) plays an essential role in el- 
ementary particle physics and quantum field theory (QFT). In QFT this phenomenon is 
well observed in Nambu - Jona-Lasinio (NJL) type models - four-dimensional models with 
four-fermionic interactions || ||. The simplest one is presented by the Lagrangian 

n G N N 

L ^ = ^^kidtp k + i^iC^^k^kf + C^^l^kf], (1) 

k=l k=l k=l 

which is invariant under continuous chiral transformations 

ipk^e^ipk ; (k = l,...,N). (2) 

(In order to apply large N-expansion technique we use here N-fermionic version of the model.) 

Since there are no closed physical systems in nature, the influence of different external 
factors on the DCSB mechanism is of great interest. In these realms, special attentions have 
been given to analysis of the vacuum structure of the NJL type models at nonzero tem- 
perature and chemical potential 0, [1[, in the presence of external (chromo-)magnetic fields 

1 Poster talk given at the " 5th International Workshop on Thermal Field Theories and Their Application" , 
10-14 August 1998, Regensburg, Germany. 
2 E-mail: kklim@mx.ihep.su 

3 This report is based on works done in collaboration with A.K. Klimenko, M.A. Vdovichenko, A.S. 
Vshivtsev and V.Ch. Zhukovskii H §. 
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|| [7], ||, with allowance for curvature and nontrivial space-time topology ||, ||. Combined 
action of external electromagnetic and gravitational fields on DCSB effect in four-fermion 
field theories were investigated in [T(| [TTJ . 



In the present paper we consider the phase structure and related oscillating effects of 
the four dimensional Nambu - Jona-Lasinio model in two cases: 1) in nonsimply connected 
space-time of the form R 3 x S 1 (space coordinate is compactified) with nonzero chemical 
potential \x and 2) in Minkowski space-time at nonzero values of /x, H, where H is the external 
magnetic field. 

1.1 NJL model at [L ^ 0. 

First of all let us prepare the basis for investigations in following sections and consider phase 
structure of the model (f ) at fi ^ in Minkowski space-time. 

Recall some well - known vacuum properties of the theory (f ) at fi — 0. The introduction 
of an auxiliary Lagrangian 

L = ipidip - ip(a! + ia 2 'j5)tp ~ ^{of + oi) ( 3 ) 

greatly facilitates the problem under consideration. (In (3) and other formulae below we 
have omitted the fermionic index k for simplicity.) On the equations of motion for auxiliary 
bosonic fields <j\p, the theory (3) is equivalent to the (1) one. 
From (3) it follows in the leading order of 1/N - expansion: 



exp(iiV S e f /(er li2 )) = J DipDip exp(i J Ld 4 x) 



where 

r 2 



/(J ~\~ o~ 
d A x 1 2 -i\ndet(id - ai - 275O2). 

Supposing that in this formula 0i j2 are not dependent from the space - time points we have 
by definition: 

S e ff(<7ip) = -^0(0X2) J d 4 x, 

where (£ = \Zof~+ cr l) : 

V (*i,2) = % + 2z f -£P l n( S 2 - p 2 ) ee 7 (E). (4) 



2G J (2tt 

Introducing in (4) Euclidean metrics (po ~^ ipo) and cutting off the range of integration 
(p 2 < A 2 ), we obtain: 

The stationary equation for the effective potential (5) has the form: 



Now one can easily see that at G < G c = An 2 /A 2 eq. (6) has no solutions apart from S = 0. 
Hence, in this case fermions are massless, and chiral invariance (2) is not broken. 

If G > G c , then Eq. (6) has one nontrivial solution S (G, A) ^ such that F(E ) = 0. 
In this case So is a point of global minimum for the potential Vo(S). This mean that 
spontaneous breaking of the symmetry (2) takes place. Moreover, fermions acquire mass 
M =£„(£, A). 

Let us now imagine that /x > and temperature T ^ 0. In this case one can find effective 
potential V^t(£) if the mesure of integration in (4) is transformed according to a rule 



Summing there over n [0] and directing in the obtained expression temperature to zero, we 
have: 

V (E) - 2 I ^9(n - v^+7)(// - V^+?), (7) 



(2tt)* 



where is the step function. Integrating in (7), we find 



16vr 2 



-2/iV^-S 2 + S 4 ln [(//+ V/x 2 -E2) 2 /E 2 ]|. 



(8) 




FIG.l Phase portrait of the NJL model at nonzero /i and for arbitrary values of fermionic mass M. Phases B and C are massive 
nonsimmetric phases, A is chirally simmetric phase. Here /12c = M, ^l\ c = u \M 2 + A 2 /M 2 ), Af2c = A/(2.21...), Mi c is 
the solution of the equation (j,f c (Mi c ) = A 2 /(4e). In the phase B particle density in the ground state is equal to zero. However, 
in phase C particle density is not zero. 
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It follows from (8) that in the case G < G c and at arbitrary values of chemical potential 
chiral symmetry (2) is not broken. However, at G > G c the model has a rich phase structure, 
which is presented in the Fig.l in terms of \x and M. (At G > G c one can use the fermionic 
mass M as an independent parameter of the theory. Three quantities G, M and A are 
connected by the Eq. (6).) In this Figure the solid and dashed lines represent the critical 
curves of the second- and first-order phase transitions, respectively. Futhermore, there are 
two tricritical points a and (3, two massive phases B and C with spontaneously broken chiral 
invariance as well as the simmetric massless phase A on the phase portrait of the NJL model 
(detailed calculations of the vacuum structure of the NJL model one can find in [|TJ). 

2. Phase structure of the NJL model at fi ^ and in the 
R 3 x S* 1 space-time 

It is well-known that unified theory of all forces (including gravitation) of the nature is 
not constructed up to now. Since in early Universe the gravity was sufficiently strong, so 
one should take it into account, many physicists study quantum field theories in space- 
times with nontrivial metric and topology. At this, NJL model is the object of special 
attention (see the review [10]]), because the idea of dynamical chiral symmetry breaking 



is the underlying concept in elementary particle physics. There is a copious literature on 
this subject || |TI| [H], |T3]]. In particular, the investigation of four-fermion theories in the 
space-time of the form R d x S 1 x • • ■ x S 1 is of great interest fl3fl . The matter is that such 



space-time topology occurs in superstring theories, in description of Casimir type effects and 
so on. 

In the present section the NJL model in the R 3 x S l space-time and at /x 7^ is considered 
since great amount of physical phemomena take place at nonzero particle density, i.e. at 
nonzero chemical potential. Here space coordinate is compactified and the circumference S 1 
has the length L. For simplisity we study only the case with periodic boundary conditions: 

ip(t,x + L,y,z) = ip(t,x,y,z). 

2.1 Phase structure 

In order to find the effective potential V^£,(£) at /j, 7^ and L 7^ 00, we need to transform 
the integration over pi in (7) into a summation over discrete values p\ n according to a rule 



n = 0,±1,±2, 



n=—oo 

The resulting expression is 

00 

V^l(S) = V L (E) --^Yl a ^ ' V^ 2 + (2vrAn) 2 )• 

ra=0 

• (/x - VS 2 + (2vrAn) 2 ) 2 (/i + 2^Y? + (2vrAn) 2 ), (9) 
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where 



oo 

V L {E) = V {E) - -4t f dxx 2 ln[ 1 - exp(-Lv / x 2 + E 2 ) ], (10) 
7r -L y 







ot n = 2 — <5„o and Vo(E) is given in (5). The stationary equation for the function (9) has the 
form 

dVyLp) dV L (E) , AE^ , 

— ^ — = g S + — 2^ a » g (^~ V s +( 2?rAn ) )• 

n=0 

2E 

•(/i - VS 2 + (27rAn) 2 ) = —0(E) = 0. (11) 

The case /i = 0, A = 1/L > 0. Putting in (9) /i equals to zero, we obtain effective potential 
V^(E) (11) in the case under consideration. This function at G > G c has a global minimum 
point Eo(A) > 0, which means that for all values of A > chiral invariance of the model is 
spontaneously broken. Obviously, Eo(A) — > M at A — > 0. 

When G < G c , E (A) = at A < A and E (A) > at A > A (here fj - ^ = ^Ag). In 
the point A = Ao there is a second order phase transition from a symmetric to a nonsymmetric 
phase of the model, because at A — > Ao+ 

Eo(A) = ^7r(A-A ) + o(A-Ao), 

i.e. the order parameter E (A) is a continuous function in the point A = Ao- At A — ► oo for 
all values of coupling constant G we have 

E (A) ~ 2ttA(2.719...). 

Details of above calculations and of following ones are presented in 0. 

The general case jtt, A ^ 0. We shall find a one-to-one correspondence between points of 
the plane (A, /x) and the phase structure of initial model. It is very convinient to divide this 
plane into regions tOk- 

oo 

fa A) = (J u k ; uj k = {(/i, A) : 2ttAA; < fi < 2n\(k + 1)}. (12) 

fc=0 

In co>o only the first term from a series in (9) is nonzero, in u\ only the first and the second 
terms are nonzero and so on. In order to obtain a phase structure one should study step by 
step the global minimum point of the function V^i(E) in regions uq,ui,.... Omiting technical 
details we show at once resulting phase portraits at G\ = (0.917...)G C < G < G c (see Fig. 2) 
and at G c < G < (1.225...)G C ee G 2 (see Fig.3) as well [|. 

One can see in Fig. 2 only two massive nonsymmetric phases B and C. In contrast, there 
are infinitly many massive phases Ck{k = 0, 1, ...) in Fig.3. In phase B particle density is 
identically zero, but in C and in all C k phases this quantity is not zero. In both figures there 
are also infinitly many simmetric massless phases A k (k = 0, 1, ..) of the NJL model. 

The line /io(0)c n C2Ci& in Fig.3 is the critical line /i c (A) of the second-order phase transi- 
tions, where chiral symmetry is restored. The /x c (A) is defined by the equation 

0(0) = 0, (13) 
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where is given in (11). Critical lines I1J2, ■■■ are solutions of equations 

0O fc ) = <t>{ y/» 2 ~ (27tA;A) 2 ) = 

for k = 1, 2, ... respectively. Boundaries between massless phases in both figures are bound- 
aries between regions u>k from (12). 

Phase structure of N JL model at another values of coupling constant G is presented in |2j , 
where one can also find more detailed description of above phase portraits at G\ < G < G 2 - 



= 27rAn 



fi = 47rA 




2ttA 




FIG. 2 Phase portrait of the _R 3 X S 1 NJL model at fi and Gi < G < G c (A = 1/L). Dashed lines are critical curves of 
first-order phase transition, solid lines correspond to second-order critical curves. Points a and b are tricritical ones. There is 
a cascade of massless symmetric phases A k (k = 0, 1, 2...). 



2.1 Effects of oscillations 

Now let us show that, due to a presence in a phase structure of the NJL model of cascades 
of massless A\~ as well as massive phases, one can observe oscillations of some physical 
parameters. We shall consider only the case G c < G < G 2 - 
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The continuous physical quantity f(x) is called as an oscillating one at x — > a, 
if there exist a monotonically encreasing (decreasing) sequence {x n } such that: 
i) x n — > a at n — > oo, ii) /(x) is a continuous function at each point x n and hi) /'(x) 
is a discontinuous function at points x n . 

At zero temperature oscillating quantity satisfies, as a rule, to this definition (see, for 

Of course, at nonzero 



example, magnetic oscillations in quantum electrodynamics |T4], |15| 
temperature one can observe a more smooth behaviour of oscillating parameters. Since we 
shall deal with zero temperature case only, the above cited definition of oscillations is well 
suited in the framework of present paper. 



Mo) 



M 




A 



FIG. 3 Phase portrait of the B? X S 1 NJL model at fi ^ and G c < G < G2 (A = l/L). Dashed lines are critical curves of 
first-order phase transition, solid lines correspond to second-order critical curves. Points a and b are tricritical points. There 
are cascades of massless symmetric phases as well as massive phases Cfe (k = 0,1,2...). The line Ma is fi = So(A) and 
ft c (P) = 2tt\ /V6, where - ^- = -^A§. 

Oscillations of the critical curve // C (A). Recall that /x c (A) is the solution of the equation 
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(13). Evidently, inside an arbitrary region ujk (see (12)) this function has the form 

27r{[6fc(fc + l) + l]A 2 + Ag} 
MA)^ = A* W (A) 6(2A;+1)A > ( 14 ) 

where Ao is given in the caption to Fig.3. Hence, 

H c (\) = //(k)(A) at tfc+i<A<t fc , fc = 1,2,3,..., (15) 

where is such a value of parameter A, that the curve /i c (A) crosses the line /i = 2nkX, i.e. 
the right boundary of uj^. 

h = ; A ° ■ (16) 

vW - 1 V ; 

Note, H(k){tk) = fi(k-i)(tk), so the function /i c (A) (15) is a continuous one at A > 0. It follows 
also from (15) that 

d£i( fc _i)(A) tt(2 - 6fc) 



rfA 



A^i fc+ 3(2/e - 1) 



< 0, 



dfi k {\) tt(2 + 6k) 



dX \-^t k _ 3(2k + 1) 

Last inequalities mean that at infinite set of points t k (k = 1, 2, ....) the function /i c (A) is not 
differentiated. According to the above given definition, the critical curve /x c (A) oscillates at 
A — > or, equivalently, at L — > oo (see Fig.3). 

Finally, let us present this oscillations of the /i c (A) in a manifest form. We need the 
following Poisson's summation formula [|T6 1 : 

oo 

oo oo „ 

^« n $(n) = 2^a fc / §(x) cos(2nkx)dx, (17) 

n=0 k=0 q 

where a n = 2 — 5„o- Using it in equation (13), one can easily find at A — > 0: 

* (A) ~ ~w 1 1 + ^a§ ^ — * — j ' (18) 

From (18) it follows that /x c (A) has an oscillating part, which oscillates at L — > oo with 
frequency A /(2a/6). 

Oscillations of the fermionic condensate. Fermionic condensate is defined as < ipip >, and 
in NJL model it equals to < £ >. Since the last quantity is the global minimum point £(//, A) 
of an effective potential, we should study nontrivial solution of the stationary equation (11). 
Detailed analysis of £(//, A) was done in p|, and this quantity at M < \i < fi c (0) and at 
A — > (L — > oo) behaves as 

A 2 ^ 2 - ™ 2 (/^) cos(ra v /// 2 - m 2 (/i) L) 2 
£(/,, A) = mfc*) + g + o(A ), (19) 



S 



where m(/i) equals to E(/i, 0), f'ijn) is the derivative of the function 



f( m ) = F(m) + ^ 2 




(20) 



and Fill) is given in (6). It is clear from (19) that S(/i, A) has oscillating part, which 
oscillates at L — > oo with frequency ^ fi 2 — m 2 (jj) / (2n). 

Oscillations of the particle density. Suppose, M < /i < /i c (0). Then the thermodynamic 
potential (TDP) A) of the NJL system is equal to the value of its effective potential 
at the global minimum point E(/z, A), i.e. fl(fi, X)= V^(E(/i, A). It is well-known that 
thermodynamic potential defines the particle density n(/i, A) through the relation: n(/i, A) = 
—dQ(fi,X)/dfi. Hence, 



n(/i, A) 



<9V^(E) «9^ L (E)<9E 



<9/i 



<9E <9^ 



E=S(/i,A) 



A 

2^ 



^ ar n e(/i - v / S 2 ( / u,A) + (27rAn) 2 )(/i 2 - S 2 (/i, A) - (27rAn) 2 ). 



(21) 



n=0 



Using in (21) the Poisson's summation formula (17) H, we see that at fi = const and L — > oo 



n(//, A) 



m 2 (/x)) 3 / 2 



3tt 2 



+ A^ 



m(/z)(/i 2 — m 2 (/z)) 



2\//i 2 — m 2 (yu) 



E 

n=l 



COS 



{n^J fi 2 — m 2 (fi)L) 



n 2 n 2 



+ o(A 2 



(22) 



where m(/i) is the fermion mass atA = l/L = (see (19)), /(tb) is defined in (20). From 
(22) one can easily see that particle density in the ground state of the NJL model oscillates 
with frequency a/ '/i 2 — m 2 (/i) / (2ir) . 

Oscillations of the pressure. Let us suppose that /i > /i c (0), i.e. for sufficiently large 
values of L the global minimum point of the effective potential equals to zero. In this case 
the TDP of the model equals to V^^(0). Then, at L —>■ oo the TDP A) oscillates with 
frequency /i/(27r), because it looks like M: 



ft(/z,A)=V £ (0)- 



127T 2 



E 

fc=0 



4A 4 



■ sm{jj,kL) — -7 cos{j2kL) 

n 2 k 



7T 



In our case the pressure in the system is defined as p = —d(LQ)/dL. Using the above 
expression for f2(/i, A), we see that the pressure in the vacuum of the NJL model also oscillates 
with frequency /i/(27r). 

One can interpret the case under consideration as the ground state of the NJL system, 
located between two parallel plates with periodic boundary conditions. The force which acts 
on each of plates is known as generalized Casimir force. Evidently, this force is propotional 
to the pressure in the ground state of the system. Hence, at nonzero chemical potential the 
Casimir force of the constrained fermionic vacuum is oscillates at L — > oo. 
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3. Phase structure of the NJL model at \i ^ and in the 
presence of external magnetic field 

In the present section we shall study the magnetic properties of the NJL vacuum. At 
/i = this problem was considered in [§, ||. It was shown in || that at G > G c the chiral 
symmetry is spontaneously broken for arbitrary values of external magnetic field H, and 
even for H = 0. At G < G c the NJL system has a symmetric vacuum at H = 0. However, 
if the external (arbitrary small) magnetic field is switched on, then for all G e (0, G c ) one 
has a spontaneous breaking of initial symmetry ||. This is a so called effect of dynamical 
chiral symmetry breaking catalysis by external magnetic field. 

The brief history of this effect is the following. First of all such property of external 
magnetic field was discovered in (2+1) - dimensional Gross- Neveu (3DGN) model |L7|, [18| . 
At H = there exist two phases in 3DGN model: one of which is a massless chirally 
invariant phase {G < G c ), and another one is a massive phase with spontaneously broken 
chiral symmetry {G > G c ). However, for each value of H ^ as well as for all values of the 
bare coupling constant G > the symmetric phase is absent in a 3DGN theory, and chiral 
symmetry is broken down \ Tf\ , |18 | Q. Of course, in |17], [18|] a special consideration was done 



for the case G < G c , where the magnetic field induces the DCSB even for weakest attractive 
interaction between fermions (the magnetic catalysis of DCSB). It turns out, that external 
chromomagnetic field is also a magnetic catalyst of DCSB [M. The influence of temperature 



and chemical potential on this effect was studied in |18|-||21||, where the restoration of chiral 
symmetry at sufficienly large values of T and \i was predicted. Later, in [22| the explanation 
of this phenomenon on the basis of dimensional reduction mechanism was found in the 
framework of 3DGN model. The magnetic catalysis takes place in four-dimensional NJL 
model H ^3], [24j] as well as in other theories, and now it is under intensive consideration 
(see, e.g. [|H], [25|, |26|] and references therein). 



(Authors of |^4[ ^(J declare that in our paper ]I3 only the "fact that external magnetic 



field enhances a fermion dynamical mass" was established. Hence, they assert that in JL| 
only the case G > G c was considered. In reality, in |I7|, [L8], [L9| the action of external 
(chromo-) magnetic field on the 3DGN model was studied for arbitrary values of bare coupling 
constant. The spontaneous breakdown of chiral symmetry was found there for all G G (0, oo), 
including the case G < G c , and even the case of arbitrary small values of G > (the magnetic 
catalysis of DCSB).) 

In the present section we continue the investigatuion of magnetic catalysis effect and this 
time turn to the consideration of the four- dimensional NJL model at H, jj, ^ 0. 

3.1 Magnetic catalysis at fj, ^ 

Let us recall some aspects of the problem at \x = 0. Using a well-known proper-time method 



27] or a momentum-space calculations [28], one can find the effective potential V#(£) of the 

4 The consideration in is performed in terms of parameter g, such that | = jj^j — ^r, where m 
is a normalization point, g(m) is renormalized coupling constant. The connection between g and G was 
established in pl|: ~ = -k — i . So, at g < and g > one has G > G c and G < G c , respectively. 
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NJL model at H ^ 0: 

E 2 eH r°° ds 
VhV) = — + — 2 J^ -exp(- S E 2 ) coth(eFs). 

After identical transformations we have 

(23) 

where 







V^E) = K (E) H 


- V H (E) 


+ Z(E), 


V (E) = 


s 2 

2G + 


1 f°° ds 

o 2 / 3 eX P 


-*£ 2 ), 




Z(E) = 


(e#) 2 
24tt 2 


/ ^?exp(-sE 2 ), 




ME) = 


Svr 2 / 


00 <is 

— exp(-sE 2 ) 
s a 


(eHs) 


coth(eifs) 



1 



2 



(24) 



The potential Vo(E) in (24) up to infinite additive constant is equal to the function (4). 
Hence, the UV-regularized expression for it looks like (5). 

The function Z(E) is also an UV-divergent one, so we need to regularize it: 

Z(E) = T (exp(- S S 2 )-exp(- S A 2 )) + ^-Ly o -exp(- S A 2 ) 

(eF)\ E 2 (eH) 2 r ds , . 2 , 

+ -ex P (-,A 2 ). (25) 



24tt 2 A 2 24tt 2 7 s 

The infinite last term in (25) contributes to the renormalization of an electric charge and 
magnetic field as well, similar as it occures in quantum electrodynamics [p7 |. 

The potential Vff(E) in (25) has no UV divergences, so it is easily calculated with the 
help of a table of integrals f29j. The final expression for Vff(E) is: 

Vh&) = Vo(E) - ^{C(-l, x) - \[x 2 - x] lnx + ^}, (26) 

where x = E 2 /(2eif), CC^^) is the generalized Riemann zeta-function and C'( — 1> X ) = 
d((u, x)/dv\ v= -\. The global minimum point of this function is the solution of the stationary 
equation: 

^(E) = ^{F&) - 1(E)} = 0, (27) 
where F(E) is given in (6), and 

7(E) = 2eH{lnT{x) - ^ ln(27r) + x - ^{2x - l)lnx} 

= / ^exp(-sE 2 )[ eifscoth(etfs) - 1] (28) 
Jo s 

For arbitrary fixed values of H, G there is only one nontrivial solution E (H) of equation 
(27), which is the global minimum point of V#(E). 
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Hence, at G < G c and H = the NJL vacuum is chirally symmetric one, but arbitrary 
small value of external magnetic field H induces the DCSB, and fermions acquire nonzero 
mass T, (H) (the effect of magnetic catalysis of DCSB). 

In the present paper we shall consider only the case G < G c . Therein, T,q(H) is a 
monotonically increasing function versus H. Apart from, at H — > oo 



and at H — > 



eH , 1 An 2 

~ — exp{ — 

7r ei/ G 



A 2 )} 



(29) 



(30) 



Now let us consider a more general case, when H ^ 0,/i ^ 0. In one of our prevoius 
papers []2T| the effective potential of a 3DGN model at nonzero H, /j and T was obtained. 
Similarly, one can find an effective potential in the NJL model at H,T,fi=£ 0: 



k=0 



dp In { [1 + exp-^+^j [1 + exp-ft*"-!*')] } , (31) 



where j3 = 1/T, a k = 2 — 5ok, e k = \/S 2 + p 2 + 2eHk, and a function Vff(S) is given in (26). 
Tending a temperature in (31) to the zero, we have the effective potential of NJL model at 

eH x 



V, 



Hfi\ 



Vh-(E) - ^ a fc / dp(p - e k )6(fi - e k ) 



(32) 



fc=0 



which can be easily transformed to the form 

eH 



fc=0 ^ 



fi+ \fli 2 - s 



(33) 



where s k = \/T? + 2eHk. Finally, let us present the stationary equation for the potential 
(33): 



d S 



F(S) - /(£) + 2e# ^ a^(/x - s fc ) In 



fc=0 



0, (34) 



In order to get a phase portrait of the model one should find a one-to-one correspondence 
between points of (/x, H) plane and global minimum points of the function (33), i.e. we need 
to solve equation (34), find a global minimum E(/i, H) for potential (33), to study properties 
of S(/i, H) versus (/i, H). 

In order to greatly simplify this problem let us divide the plane (/i, H) into a set of regions 

u k : 



(p, H) = |J f; fc = {(p, If) : 2eFA; < /i 2 < 2eH(k + 1)}. 



(35) 



fc=0 



In the region ujq only a first term is nonzero from a series in (34-35). So, one can find that 
for the points (fi,H) G uj which are above the line L={(/i, if) : \i = T> (H)}, the global 
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minimum is at the point £ = 0. Just under the curve L the point £ = Y>o(H) is a local 
minimum of the potential (33), and £ = H (H) transforms to the global minimum when 
(fi,H) lies under the critical curve of the first order phase transitions \i = /i c (if), which is 
defined by the following equation: 

VhM = V Hli (E (H)). (36) 
In the region u one can easily solve this equation: 

2-7T 

Hc{H) = -j=Wh(0) ~ V H {MH))] 1/2 - (37) 
yen 

Hence, we have shown that at /i > ^ C {H) (G < G c ) there is a massless symmetric phase 
of the NJL model (numerical investigations of (33-34) gives us the zero global minimum 
point for the potential Vj/^E) in other regions ui,U2, ■■■ as well). The external magnetic 
field ceases to induce the DCSB at fj, > fJ, c (H) (or at sufficiently small values of magnetic 
field H < H c (fj,), where H c (/j,) is the inverse function to /i c (H)). But, under the critical curve 
(37) (or at H > H c (fi)) due to a presence of external magnetic field the chiral symmetry is 
spontaneously broken. Here magnetic field induces dynamical fermion mass £ (i?), which 
is not /i-dependent value. 

At last, we should remark that in the NJL model the magnetic catalysis effect takes 
place only in the phase with zero particle density, i.e. at /i < fi c (H). If /i > /i c (H) we have 
symmetric phase with nonzero particle density, but here the magnetic field can not induces 
DCSB. 



3.2 Magnetic oscillations 

In a previous case we have shown that points (fi,H), which are above critical curve fi = 
fi c (H), correspond to the chirally symmetric ground state of the NJL model. One fermionic 
excitations of this vacuum have zero masses. At first sight, properties of this simmetric 
vacuum are slightly varied, when parameters \i and H are changed. However, it is not 
so and in the region \x > fi c {H) we have infinitly many massless symmetric phases of the 
theory as well as a variety of critical curves of the second order phase transitions. On the 
experiment this cascade of phases is identified with oscillations of such physical quantities 
as magnetization and particle density. Let us prove it. 

It is well-known that a state of the thermodynamic equilibrium (= the ground state) of 
arbitrary quantum system is described by the thermodynamic potential (TDP) Q, which is a 
value of the effective potential in its global minimum point. In the case under consideration 
the TDP H) at fj, > fJ, c (H) has the form 

Q(fi,H) = V HlI (0) = V H (0)- 

oo 

- Yl akd ^ ~ ^\l^ 2 - e l - e l ln[(^ 2 -e| + fi)/e k ] }, (38) 

fc=0 

where = y2eHk. We shall use the following criterion of the phase transitions: if at least 
one first (second) partial derivative of f2(/i, H) is a discontinuous function at some point, 
then it is a point of the first (second) order phase transition. 
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Using this criterion let us show that boundaries of uo^ regions (35), i.e. lines 4 = {(/x, H) : 
\x = y2eHk} (k = 1,2,...), are critical lines of second order phase transitions. In arbitrary 
region uj^ the TDP (38) has the form: 



8=0 



From (39) one can easily find 



k-l 



dfx 



as well as: 



d 2 Qi 



d 2 n 



fc-i 



{n,H)-*l k ^ 



2<k 2 JW^\ 



-oo. 



(39) 



(40) 



(41) 



The equality (40) means that the first derivative dQ/d/j, is a continuous function at all lines 
Ik. However, the second derivative d 2 fl/(dfj,) 2 has an infinite jump at each line 4 (see (41)), 
so these lines are critical curves of second order phase transitions. (Similarly, one can prove 
the discontinuity of d 2 VL/(dH) 2 and d 2 £l/dfidH at all lines l n .) 

Let the chemical potential be fixed, i.e. \x = const. Then on the plane (/i, H) we have 
a line, which crosses critical lines ■■■ at points H\,H2, ... correspondingly. The particle 
density n and the magnetization m of any thermodynamic system are defined by the TDP in 
a following way: n = —dQ/dfi, m = —dVt/dH. At fi — const these quantities are continuous 
functions over external magnetic field only, i.e. n = n(H), m = m(H). We know that all 
second derivatives of H) are discontinuous at every critical line l n . So, functions n(H) 
and m(H), continuous at the interval H e (0, oo), have derivatives broken at infinite set 
of points Hi, Hk, .... According to the definition given in section 2.1 particle density and 
magnetization oscillate at H — > 0. 

In order to present oscillating parts of n(H) and m(H) in a manifest analitic form, we 
shall use the technique elaborated in []nj, where a manifest analitic expression was found for 
oscillating part of fi(/x, H) for free relativistic electron - positron gas. This technique may 
be used without any difficulties in our case as well. 

Hence, one can rewrite the TDP (38) in the following form: 



Q((JL, H) = Qmon(fJ>, H) + Q osc (fi, H), 



(42) 



where (u = fj, 2 /(eH)): 



V mon = Vff(0) 



fi 4 {cHf 



12tt 2 4tt 3 



OO j 



(43) 



k=l 



a. 



4^3/2 



°° feH\ 3/2 

( — — j lQ(nkv) cos(7rA;z/ + vr/4) + P{ixkv) cos(nku — 7r/4)]. 
k=i ^ 71 ' 



(44) 
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(To find (43-44) it is sufficient to tend the electronic mass to zero in the formula (19) from 
flT5fl .) Functions P(x) and Q(x) in (43-44) are connected with Fresnel's integrals C(x) S(x) 

m\- 



C(x) 



S(x) 




+ \ — [P(x) sinx + Q(x) cosx] 
— 4 / — [P(x) cosx — Q(x) sinx] 



They have at x — > oo following asimptotics [30| : 



P(x) = x 



3 

--x~ 3 + ..., Q{x) 



1 -2 15 _ 4 

= — x 2 H x 4 + ... 

2 8 



The formula (44) presents the exact oscillating part of the TDP (38) for the NJL model at 
G < G c . Since in the present case the TDP is propotional to the pressure of the system, one 
can conclude that pressure in the NJL model oscillates, when H — > 0. It follows from (44) 
that frequency of oscillations at large values of a parameter (eif) -1 equals to /x 2 /2. Then, 
starting from (44) one can easily find manifest expression for oscillating parts of n(H) and 
m(H). These quantities oscillates at H — > with the same frequency /x 2 /2 and have a rather 
involved form, so we do not present it here. 

The author is grateful to Prof. U. Heinz and to the Organizing Committee of the "5th 
International Workshop on TFT and Their Application" for kind invitation to attend this 
Workshop as well as to the Deutsche Forschungsgemeinschaft for financial support. This 
work was supported in part by the Russian Fund for Fundamental Research, project 98-02- 
16690. 
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